DECOHERENCE IN A QUANTUM WALK ON THE LINE 
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- - - We have studied how decoherence affects a quantum walk on the line. As expected, 

it is highly sensitive, consisting as it does of an extremely delocalized particle. We 

' obtain an expression for the rate at which the standard deviation falls from the 

^ " ^ quantum value as decoherence increases and show that it is proportional to the 

, number of decoherence "events" occuring during the walk. 

> ■ 

; 

, 1 Introduction 

OO 

\ Quantum walks on discrete lattices are of interest to the quantum information 

I ■ community because classical random walks underpin important classical al- 

^ I gorithms. The aim is to find more powerful quantum counterparts. First it is 

' useful to understand the basic properties of quantum walks on simple lattices, 

' Two versions have been studied, continuous timeB and discrete-time coineda 
quantum walks. Here we consider only the coined version. Coined 
tum walks on a line, cycle, and the hypercube have been solved exact! 

p\J I and numerical studies used to explore quantum walks on two dimensional 

I latticedaiD. The quantum walks show a speed up over the equivalent classical 

(-H , random walk, a quadratic increase in standard deviation for a walk on a line, 

^ il a quadratically faster mixing time of a walk on a cycle and an exponentijal 

1 . increase in the hitting time to the opposite corner for a walk on a hypercuba3. 
C ■ 

^ : 

^ 2 Qauntum walk on a line 

FT' 

^ ' A quantum coin is a d-dimensional quantum system which is "flipped" by 

applying a unitary operator C. For a quantum walk on a line, d — 2 ("left" 
and "right" ) and without loss of generality C can be taken to be a Hadamard 
{tt/2) rotation. The particle then moves in a superposition of the two possible 
directions indicated by the coin, S|a;,a) = |a; + a,a), where x is the particle 
position and a e is the state of the coin. The unitary operator 

for each step of the walk is thus U = S • (C (K) I). The moments have been 
calculatedEl, for asymptotically large times T for a walk starting at the origin, 



(1 - 1/V2)T' = a'{T), (1) 
independent of initial coin state, and 

{x)a - a{l - l/V2)T, (2) 
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where a E {1, —1} is the initial coin state. The standard deviation (from the 
origin) a{T) is thus hncar in T, in contrast to Vt for the classical walk. 

3 Decoherence in the walk on a line 

To model decoherence in this system we use a discrete master equation 

pit+1) = {l-p)Up{t)U^ +pY,JP^Up{t)U^IP^, (3) 

i 

where the summation runs over the dimensions of the Hilbcrt space on which 
the decoherence occurs, either the coin, the particle, or both. The projectors 
IPi act in the computational basis, and p is the probability of a decohecence 
event happening per time step. This equation was evolved numericallytl for 
various choices of IP. All choices produced the same general form for the 
decay of a{T,p) from the quantum [p = 0) to the classical {p = 1) value, 
with small differences in the rates. The slope of a{T,p) is finite as, p ^ Q and 
zero aX p = 1. Here we calculate a{T,p) analytically for pT ^ 1 and T ^ 1 
for the case where IP is the projector onto the preferred basis {|a,x)}, i.e., 
decoherence affecting both particle and coin. 

The probability distribution for finding the particle in the state \a,x) in 
the presence of decoherence can be written, 

P(x, a, T, p) = {l- pfP{x, a, T) + p{l - pf-'P<^''> (a:, a, T) + . . . , (4) 

where P{x, a, T) is the distribution obtained for a perfect walk and 
P^'^\x,a,T) is the sum of all the ways to have exactly i noise events, e.g., 

T 

P(i) {x, «, r) = E E E P^y^ ^' ^^^Vi^^ a,T-t), (5) 

t=l y b 

where Pyb{x, a,T — t) is the distribution obtained by a perfect (no noise) 
walk starting in state \y,b) for T — t steps. For the ideal walk, cr^iT) = 
Tlx Ta x^P{x, flj T), and for the walk with decoherence, 

a2(r,p) = ^^a:2p(a;,a,T,p). (6) 

X a 

Taking Eqn. (jj) to first order in p, and substituting along with Eq. (|^) into 
Eq. (i) gives 

a^{T,p) ~ ^ I (1 - pT)Pix, a, T) + ^ F(y, b, t)Pyk{x, a,T~t)\. 

x,a I t—1 y,b J 

(7) 

The first term on the r.h.s. is (by definition) (1 — pT)a'^{T). Noting that 
Pyh{x, a, T—t) is a translation of a walk starting at the origin, Pyh{x, a, T~t) = 
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Pob{x — y,a,T — t). Relabelling the summed variable x to (x + y) then enables 
the sums over x and a to be performed in the second term, 



E E ^(2^' ^) E E(^ + y)'Pobi^^ ^^T^t) 

t—l y b X a 

T 

= P E E E P(y^ ^' ^) {''obiT -t) + 2y{x)g''^ +y^]. (8) 
t=i y b 

From Eqn. (|^), (Tq^(T — t) does not depend on b, so the summation over y 
and 6 may be performed trivially. The summation applied to y^ gives crg(i) 
by definition. This leaves only the evaluation of 

T T 

2pE E P(y^ t)y{x)lg-'^ = 2p(l - 1/V2) ^(T - t) J] y& P{y, b, t), (9) 

t=l y,b t=l y,b 

where we have used Eqn. (^. We note that this term does not depend on 
whether the initial coin state is plus or minus one and so we may include both 
these possibilities equally. Also, by the symmetry of the walk, it is possible 
to rewrite a probability function for travelling from state |0, c) to b) in the 
reverse order, i.e. as a probability for moving from 6) to |0, c). Care must 
be taken to ensure that the signs of each term due to the coefficient yb in the 
summation are maintained. We obtain 

T 



2pY^Piy,b,t)y{:. 



/Ob 



P{1 - \I^)Y,{T - t)Y, (1 - 26b,c)yh PybiO, c, t) 



t.y,b 



y,b,c 



= p(i-i/V2)^(r-<) 



^ yb Py^biO, c,t)~2j2yb PybiO, b, t)] 

y.b,c y,b 



(10) 



treating the two parts with and without a delta function independently. Ex- 
panding the summations over b and translating the particle basis by — y gives 

T 

2pEE^(y'^'*)y(^)r"'^ 



/Ob 



t=l y,b 

T 

= p{l-l/V2)J2{T-t) 



Y,yPo-i{y^c,t)-Aj2y Py,+i{y,+ht) 



(11) 



The final summation over y may be bounded above by noting that 

Fo-i(y,-i,t) = l(2/,-i|f/*|0,-i)p 

= + 1, l|[/*-i|0, -1) - (y + 1, -l|f/*-i|0, 



< i^Po,-i(y + i,c,i-i). 



(12) 



3 




probability of noise event (p) 



Figure 1. Standard deviation cr(T,p) of the particle position on the line for T = 20-500. 
Inset shows small p region scaled by converging on value of 0.06378 .... 



Using this in Eqn. (^Tj) gives, 

T 



t,y,b 



i = l 
T 



<p{l-l/V2)Y,iT~t) {1-1/V2)t + V2 ,(13) 



where Eqn. has been used for the average values. Combining these resuhs 
in the fuU expression for a^{T,p), Eqns. (R, and performing the summations 
over t using J2 * = ^iT + l)/2 and J2 1 = T^/^ + T^f^ + T/6, gives 



a\T,p)<aliT) 



l^^pT + p{V2-l) 
6 



(14) 



Taking the square root gives as an upper bound on the standard deviation. 



a{T,p) < a{T) 



1 



6V2 V2 



^(l-l/V2) + 0(p^l/T) 



(15) 



4 Summary 



Equation (|T^) compares well with simulation data, once a second order cor- 
rection for cr(T) = (1 - l/V2)^/2(r - l/T) is taken into account. The l/T 
form of this correction was found in Ambainis ct. ala, however, we determined 
the coefficient (1 — l/\/2)^/^ numerically for a wide range of evolution times 
T. The bounding procedure applied here is reasonably accurate, numerical 
studies give the coefficent of p in the above expansion as 0.09566, compared 
with the bound of 0.20711. The first order dependence is thus proportional to 
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pT, the number of decoherence events during the whole quantum walk. For 
a given decoherence rate p, the standard deviation initially decreases linearly 
in T.|_Elelated results for decoherence on the coin only can be found in Brun 
et. allia. 
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